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Abstract 

Lattice modulation and magnetic structures in magnetoelectric compounds Tbi-^Dy^MnOa 
have been studied around the ferroelectric (FE) Curie temperature Tq by x-ray and neutron diffrac- 
tion. Temperature-independent modulation vectors through Tq are observed for the compounds 
with 0.50 < x < 0.68. This indicates that ferroelectricity with a polarization (P) along the c axis 
in the i?Mn03 series cannot be ascribed to such an incommensurate-commensurate transition of an 
antiferromagnetic order as was previously anticipated. Neutron diffraction study of a single crystal 
with x = 0.59 shows that the FE transition is accompanied by the transformation of the Mn-spin 
alignment from sinusoidal (collinear) antiferromagnetism into a transverse spiral structure. The 
observed spiral structure below Tq is expected to produce P along the c axis with the 'inverse' 
Dzialoshinski-Moriya interaction, which is consistent with the observation. 

PACS numbers: 75.25,+z, 61.10.Nz, 75.80.+q 
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Research on coupling between ferroelectric (FE) and magnetic orders has revived after 
a large magnetoelectric (ME) response was reported in several compounds. In particular, 
gigantic ME effects in rare-earth manganites i?MnC>3 (R =Gd, Tb, or Dy) with orthorhom- 
bically distorted perovskite structure have stimulated considerable interest since Kimura et 
al. reported that an electric polarization P can be flopped by applying a magnetic field 
in TbMnOs.yi Evidence has been increasing that this unique ME effect is inherently con- 
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nected to the Mn-spin ordering structure. |2J, La U |5|] TbMnOs and DyMnC>3 undergo an 
antiferromagnetic (AFM) transition with a temperature-dependent modulation vector (0 
<?Mn 1) around 40 K. Hereafter, we use the Pbnm orthorhombic setting. With further cool- 
ing, the gMn value appears to be locked, and simultaneously, FE polarization appears along 
the c axis. GdMnC>3, whose orthorhombic distortion is less than that of TbMnOs, shows a 
magnetic-field-induced FE state with P || a. (5, If], Q] Synchrotron x-ray diffraction measure- 
ments of GdMnO"3 and TbMnO"3 in magnetic fields revealed that the modulation wavevector 
changes with the magnetic-field-induced electric transitions. [s| 

There are two possible mechanisms of spin-order-driven ferroelectricity (see Fig.QJ. It was 
suggested that the ferroelectricity in i?MnO"3 can be attributed to the symmetric exchange 
striction at the earlier stage, due to the locking of gMn at the FE Curie temperature Tq. 
The superexchange interaction via an anion X between two magnetic cations (M) can be 
modified with an M-X-M angle. P| Conversely, an M-X-M bond angle should be more or 
less affected by the spin arrangement. In perovskite rare-earth manganese oxides i?Mn03, 
Mn-O-Mn bonds are not straight due to rotation and tilt of the Mn06 octahedra. Therefore, 
the Mn-O-Mn-O- bonding along (110) forms a zigzag chain, as schematically shown in 
Fig. I^a). Below the AFM transition temperature Tjv, each Mn-O-Mn angle may be further 
influenced by the alignment of Mn moments. The modification of the oxygen positions with 
an AFM spin arrangement is demonstrated in an exaggerated manner in Fig. IHb), where 
the oxygen ions are represented first by dotted circles and then by solid ones. Ferromagnetic 
superexchange coupling between two neighboring Mn ions in the basal plane becomes weaker 
as the Mn-O-Mn bond angle decreases. 4 , 0] The inverse interaction should decrease the 
angle of the Mn-O-Mn bonds between antiferromagnetically aligned neighboring Mn spins, 
thereby resulting in a lattice modulation with a propagation vector = 2gMn- in the original 
zigzag chain, the oxygen ions between up and down spins move in the same direction when 
the magnetic modulation vector is represented as 1/n, where n is an even integer. More 
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generally, a periodic modulation of Si ■ Sj with a modulation vector of m/n can induce 
ferroelectricity in the chain, when m and n are odd and even integers, respectively. 

As an alternative mechanism of the spin-order-driven ferroelectricity, it has recently been 
pointed out that transverse spiral spin ordering can induce ferroelectricity regardless of 
the commensurability between spin modulation and the lattice, as shown in Fig. HTc). (ill 



111 ] This mechanism is closely related to the antisymmetric superexchange or the so-called 



Dzialoshinski-Moriya (DM) interaction. |12| A bending Mi-X-Mj bond lacking an inversion 
center favors a canted spin arrangement. This antisymmetric superexchange term is linear 
to the spin-orbit coupling term, and is expressed as D • (Si x Sj), where D is a constant 
vector. By considering the inverse effect, the stable position of the X ion between two 
magnetic M ions can be modified using the vector product of the two magnetic moments, 
Si x Sj. Since the vector product of any two neighboring M moments is constant in case 
of transverse spiral ordering, all the X ions are displaced in the same direction. 

Recent neutron diffraction studies have demonstrated a change in the magnetic structure 
of TbMn0 3 at T C .Q,Q Kenzelmann et al. discussed the magnetic symmetry by means 
of a model fitting of lq components of magnetic reflections, which indicated a sinusoidal 
(collinear) and spiral (non-collinear) spin ordering above and below Tc, respectively. From 
their precise measurements, it was also found that the modulation vector of antiferromag- 
netism is nearly locked at temperature greater than Tc by a few K.^J These results are 
inconsistent with the exchange striction mechanism. However, the spiral spin structure could 
not be evidenced clearly because the higher Fourier component was neglected. 

This paper reports that it is the spiral spin ordering and not the locking of gMn that 
plays a dominant role in producing ferroelectricity in RMn0 3 . Up to T/v and even above 
Tc, solid solution of TbMnC>3 and DyMnC>3 in an appropriate ratio shows long- wavelength 
antiferromagnetism with a temperature-independent modulation vector. A neutron diffrac- 
tion study on Tbo.4iDy .5gMn03 with c/Mn ~ 1/3 enables us to unambiguously determine 
the spin structure without a model fitting. The results clearly indicate a spiral magnetic 



<?Mn (~ 1/3). 

4] Off-resonant 



ordering below Tc, which changes from a collinear structure with an identica 

Single crystals of Tbi-^Dy^MnOs were grown by a floating zone method, 
single- crystal x-ray diffraction measurements were performed on Beamline 4C at Photon 
Factory, JAPAN. A 13-keV monochromatic x ray was focused on a single crystal in a closed- 
cycle He refrigerator, which was mounted on a four-circle diffractometer. 



3 



(a) 



a b 



Mn 




(b) 



4(a+b) 



FIG. 1: (Color online) Two possible mechanisms of spin-order-driven ferroelectricity in orthorhom- 
bic i?Mn03. (a) A zigzag Mn-0 chain along (110), i.e., in the basal plane (001). All the bonds 
bend with the GdFeO"3-type distortion above the Neel temperature T/v- Small and large balls 
represent Mn and O ions, respectively, (b) Ferroelectricity induced by antiferromagnetic ordering 
of Ising-type Mn spins with q = 1/4 and the resultant exchange striction. Oxygen ions located 
between up and down spins are displaced in the same direction. A similar situation can be expected 
when the modulation vector of the Mn-spin moments (not necessarily collinear) is m/n, with odd 
m and even n. (c) All the oxygen ions are displaced in the same direction with a transverse spiral 
alignment of the Mn-spin moments. 

First, we show that the lattice modulation vector q^ can be locked even above Tc in 
some solid solutions of TbMn03 and DyMnOa. Figures Efa) and (b) show the evolution of a 
superlattice x-ray diffraction peak with temperature in Tbi-^Dy^MnOs with x = 0.50 and 
0.68. A superlattice reflection is observed at (0 2.66 3) and (0 3.32 1) for x = 0.50 and 0.68, 
respectively, and the qi value can be regarded as 0.66 and 0.68, respectively. The position 
of the superlattice reflection is not dependent on the temperature for both the compounds. 
In Figure Efc), the qi values for some compounds are plotted against temperature. The Tc 
values are also shown by arrows in the panel. Here, it should be noted that the q^ value 
should be equal to twice the magnetic modulation vector of the Mn ions, <7Mn-Q| Therefore, 
we can conclude that the nearly commensurate qun values (~ 1/3) of the compounds with 
x = 0.50 and 0.68 are independent of temperatures below T/v within the resolution limit, 
in contrast to TbMnC^ and DyMnCV This result contradicts the simple scenario in which 
the commensurate locking of the AFM wavevector causes an FE transition in orthorhombic 
i?Mn0 3 . 
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FIG. 2: (Color online) (a) and (b) X-ray diffraction profiles of (0 k 1) scan in Tbi-^Dy^MnOs 
with x = 0.50 and 0.68 at various temperatures. A peak at k ~ 2.63 for x = 0.50 appears 
even above Tjv and can be ascribed to x-ray diffraction from a Beryllium window, (c) Lattice 
modulation wavenumber qi in some Tbi_ x Dy x Mn03 compounds as a function of temperature. 
Here, qi = 2qy^ n (see text). Arrows indicate the ferroelectric Curie temperatures. 

Another possible origin of ferroelectricity in the orthorhombic i?Mn03 system is a change 
in the spin structure as aforementioned. Hence, the determination of the magnetic structure 
in the Tbi.^Dy^MnOs compound provides us important information for discussing a mech- 
anism of ferroelectricity in the system. It is very difficult to carry out a model-free analysis 
of a magnetic structure with an incommensurate wavenumber because there are infinite in- 
dependent magnetic sites in an incommensurate antiferromagnet. Therefore, we selected a 
compound with x = 0.59, where the value can be regarded as 1/3, judging from the gMn 
values corresponding to x = 0.50 and 0.68, as shown in Fig. |21 The measurements of the 





5 



position of the magnetic peaks showed that the gMn value of the compound with x = 0.59 
was also independent of temperature (not shown in with a figure) within the resolution limit 
of the present neutron diffraction measurement. T/v and Tq of the compound are 41 K and 
22 K, respectively. The AFM ordering of rare-earth moments occurs at 7 K. 

Single-crystal neutron diffraction measurements were carried out on a four-circle neu- 



tron diffractometer (FONDER) instal 
Energy Research Institute, Japan 



ed in the guidehall of JRR-3M at the JAPAN Atomic 
151 ] The FONDER diffractometer was operated with a 
wavelength A = 1.24 A that was obtained by using a Ge (311) monochromator. A single 
crystal with a diameter of ~ 2.5 mm and a length of ~ 7 mm was used for the neutron 
measurements. Integrated intensities were collected for 176 nuclear reflections at 15 K and 
150 and 67 magnetic reflections at 15 K and 30 K, respectively. (At 30 K, we did not 
perform measurements for high-angle magnetic reflections where no peak was discernible at 
15 K.) The absorption was corrected by using DABEX software. The atomic coordination 
and isotropic thermal parameter of each atom as well as the scaling factor and extinction 
were determined from the intensities of the nuclear reflections by using a least-square pro- 
gram. Here, we neglected the superlattice modulation of the atomic positions because it 
was too weak to be detected by neutron diffraction. The magnetic moments were calculated 
from the intensity data of the magnetic reflections, including the magnetic contribution at 
fundamental reflections, by using a Levenberg-Marquardt-type least-square method. The 
magnetic form factors of Mn 3+ , Tb 3+ , and Dy 3+ were assumed to be isotropic. 3] At 30 
K, the fitting parameters were three components at each of the twelve Mn (m^, m^, and 
mj. c ). At 15 K, the magnetic moments at rare-earth sites were also considered to fit the 
data. Many trials of gradient descent algorithms with different initial guesses verified that 
the obtained result was one of the best minima. Although there are several solutions with 
almost similar residuals, they are essentially identical: they can transform into each other 
with some symmetry operation. Figure 01 presents the comparison between the observed 
(^obs) and calculated (J ca i) intensities of the magnetic reflections. The reliability factors 
based on integrated intensities, Rj, are 0.092 and 0.034 for 15 K and 30 K, respectively. The 
average deviation, defined as (1/N)J2 (7 ca j — I bs) 2 /c?> where N denotes the data number, 
is satisfactorily small for both temperatures (0.7 and 1.1 for 15 K and 30 K, respectively). 

The obtained magnetic structures at 15 K (< Tq) and 30 K (> Tq) are shown in Fig. 0] 
The a component of any Mn spin is smaller than 1 fis at both the temperatures. The b 
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FIG. 3: (Color online) Comparison between observed and calculated intensities of magnetic re- 
flections at (a) 15 K and (b) 30 K. The calculation was based on the magnetic structures shown 
in Fig. 4. 

components of the Mn spins are sinusoidally modulated along the b axis with an amplitude 
of 3.7 and 2.7 /xb at 15 K and 30 K, respectively. The clearest difference is observed in the 
modulation of the c components. The amplitude of sinusoidal modulation is 1.4 ~ 1.6 /xb at 
15 K, while it is less than 0.4 /xb at 30 K. Below Tq, the phase of the modulation of the c 
components is shifted by 0.457T-0.477T from that of the b components. Note that the phase 
shift is ±7r/2 in the case of an ideal spiral structure. This transverse-spiral alignment of spins 
breaks mirror reflection symmetry normal to the c plane, in accordance with the P || c state 
below Tq- Further, it should be noted that the trajectory of the Mn spins is not a circle but 
an ellipsoid, as shown in the inset of Fig. EJc). The value of Si x Si+i and the displacement 
of oxygen ions due to the inverse DM interaction should also be modified with 2qu n - This 
is consistent with the observed relation = 2gMn- In contrast to the spin structure at 15 
K, the difference in phase of the spin modulation between the b and c components becomes 
considerably smaller at 30 K. The calculated difference in phase is +0.207T and — 0.177T for 
sites with z = and z = 0.5, respectively. The spin moments for Mn and Tb/Dy above 
and below Tq projected onto the (100) plane are schematically drawn in Figs. HJc) and 
(d), respectively. The induced moments at the rare-earth sites are approximately 1 /xb at 
15 K. These figures clearly indicate that the FE transition is accompanied by a spin-state 
transformation from the sinusoidal collinear to transverse spiral type. 

In conclusion, synchrotron x-ray diffraction studies of the (Tb,Dy)Mn03 system ex- 
cludes the scenario in which an incommensurate-to-commensurate locking of the modu- 
lation wavevector of an antiferromagnetic Mn-spin ordering would induce ferroelectricity 
with P || c. A model-free analysis of the spin structure of system with gMn ~ 1/3 indicates 
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FIG. 4: (Color online) (a) and (b) Alignment of Mn-spin moments at 15 K and 30 K obtained 
from least square fitting of the intensity of the magnetic reflections. Note that there are two basal 
Mn-0 planes in a unit cell with z% = and 0.5. Lines show the best fits with sinusoidal functions, 
(c) and (d) Schematic drawing of the magnetic moments projected onto the (100) plane at 15 K 
and 30 K. Right panels show the trajectories of the Mn spins. 

that the ferroelectric transition results from the sinusoidal-to-spiral spin-structure transi- 
tion. The close relation between non-collinear spin ordering and electric polarization will 
result in new methods for discovering new magnetoelectric materials. 
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